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Abstract 

We measure the extinction law in a galaxy's spiral arm and interarm regions using a 
visual and infrared [BVRJHR) imaging study of the interacting galaxies NGC 2207 and 
IC 2163. This is an overlapping spiral galaxy pair in which NGC 2207 partially occults 
IC 2163. This geometry enables us to directly measure the extinction of light from the 
background galaxy as it passes through the disk of the foreground galaxy. We measure the 
extinction as a function of wavelength, and find that there is less extinction in the optical 
bands than expected from a normal Galactic extinction law. This deviation is significantly 
larger in the interarm region than in the spiral arm. The extinction curve in the spiral 
arm resembles a Milky Way Ry = 5.0 dust model and the interarm extinction curve is 
flatter ( "greyer" ) still. We examine the effect of scattering of background galaxy light into 
the line of sight and find that it is negligible. We also examine the effect of an unresolved 
patchy dust distribution using a simple two-component dust model as well as the clumpy 
dust model of [Witt fc Gordon (1996)| . Both models clearly demonstrate that an unresolved 



patchy dust distribution can flatten the extinction curve significantly. When fit to the data, 
both models suggest that the observed difference between the arm and interarm extinction 
curves is caused by the interarm region of NGC 2207 having a higher degree of dust 
patchiness (density ratio between high-density and low-density phases) than the spiral arm 
region. We note that an unresolved patchy dust distribution will cause us to underestimate 
the average column depth of gas in a galaxy if based solely on the visual extinction. It is 
much better to use the infrared extinction for this purpose. 



Subject headings: galaxies: dust, extinction, scattering — galaxies: individual 
(NGC 2207, IC 2163) 
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1. INTRODUCTION 



Little is known about the extinction law for the large-scale dust distribution in spiral 
galaxies. There is even some controversy as to how large the total extinction through a 
spiral galaxy is. [Valentiin (1990)| presented evidence that spiral disks are opaque. Purstein 



let al. (1991)1 pointed out selection effects in Valentijn's study; nevertheless, they also found 
that galaxies are optically thick. On the other hand, these results are in direct conflict with 
the simple fact that we can see out of the disk of our galaxy. 

A more direct approach to solving this problem has been taken by [White fc Keel| 
(1992)] . They proposed that the optical depth in a spiral galaxy may be measured directly if 



it lies in front of another object. The ideal case consists of a face on symmetric spiral galaxy 
half overlapping with a symmetric, background spiral or elliptical galaxy. The galaxies need 
to be symmetric so that their non-overlapping parts may be used to estimate the surface 
brightness of their overlapping parts. [White fc Keel (1992)| obtained B and / images of such 



a candidate, AM1316-241, and found that the observed extinction is largely confined to the 
spiral arm regions, leaving the interarm regions mostly transparent. Most recently, [White 



[et al. (1996)[ have studied 10 overlapping galaxy pairs and reached a similar conclusion in 



all cases. The advantages of this direct approach to studying the opacity of spiral disks 
are numerous and are described in detail by [White et al. (199"B}]. Because the background 



galaxy is behind the foreground galaxy, the approximation of the foreground galaxy as a 
dust screen is a good one. This geometry is much simpler than that in a normal disk galaxy 
where dust and stars are mixed. 

Previous studies of occulting galaxies have been restricted to visual wavelengths. In 
this paper we perform a similar study of the extinction through the disk of a spiral galaxy, 
extending our wavelength base into the near infrared. We present a six band [BVRJHK] 
study of the interacting galaxy pair NGC 2207 and IC 2163. This pair consists of two 
symmetric, almost face-on interacting spiral galaxies which are partially overlapping. The 
galaxies' Hubble types are SAB(rs)bc pec and SB(rs)c pec, for NGC 2207 and IC 2163, 
respectively ( |de Vaucouleurs et al. 199"l| ). These galaxies have been observed and modeled 



in detail by D. Elmegreen et al. (1995) and B. Elmegreen et al. (1995] , whose HI and 



optical observations indicate that the two galaxies have probably gone through a recent 
nonmerging encounter, p. Elmegreen et al. (1995][ do not find any significant excess of star 



formation in these galaxies as a result of their interaction. 

By determining the extinction through the foreground galaxy disk at many wavelengths, 
in both spiral arm and interarm regions, we can sample the effective large-scale extinction 
curve, and compare it to Galactic dust models. In sec. 3 we describe how we derive the 
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extinction curve. In sec. 4 we present our results and discuss deviations of our observed 
extinction curve from Galactic dust models. 



2. OBSERVATIONS 



NGC 2207 and IC 2163 are good candidate galaxies for the following reasons: 1) both 
galaxies display good symmetry, without which we cannot estimate their intrinsic colors in 
the overlapping region; 2) The foreground galaxy is close to face on; this geometry provides 
a small optical path through the foreground disk and facilitates our estimation of scattering 
effects since the optical properties of the disk are dependent on inclination with increasing 
sensitivity at high inclinations ( Pruzual et al., 1988|) ; 3) the two galaxies are positioned so 
that they partially overlap, but still have large enough non-overlapping regions to allow 
estimates of their true colors; 4) The simulations of p. Elmegreen et al. (199"5]| supply us 



with a distance between the two galaxies which is necessary for determining the effect of 
scattered light on our extinction estimates. 

These galaxies were observed as part of an imaging survey being carried out by the 
Ohio State University of ~ 200 galaxies. The survey's goal is to produce a library of 
photometrically calibrated images of late-type galaxies from 0.4 to 2.2/im. For notes on the 
observation and reduction techniques see Progel et al. (1997)| , or for individual examples. 



see IQuillen et al. (1994)| , and [Quillen et al. (1995)| . All the images were obtained at CTIO 



on the 1.5m telescope. The BVR images were observed on 1995 March 7 using the Tek#2 
10242 pixel CCD with a spatial scale of 0."44/pixel. Total on-source exposure times were 
20min, 15min and lOmin for B, V, and R respectively. The JHK images were observed 
on 1995 March 14 using the NICMOS 3 256^ pixel infrared array with a spatial scale of 
l."16/pixel. Total on-source exposure times were 16min, 15min and 29min for J,H, and K 
respectively. 

A color image was constructed and is displayed as a color plate in Figure 1. The 
color image was composed by combining the B,R, and if-band images, coding them as 
Blue, Green, and Red intensities, respectively. The relative scaling between the images was 
chosen so as to neutralize the colors of most field stars, which serves to emphasize dusty 
regions as reddish and regions of star formation as blue, while the intensities preserve the 
general photometric appearance of the galaxies. 

It is immediately noticeable (see Fig. 1) that the background spiral galaxy can be 
easily seen through the foreground galaxy even though a large portion of it lies directly 
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behind one of the foreground galaxy's spiral arms. Moreover, reddening of the background 
galaxy is plainly visible in this region, and the dust content of the foreground galaxy, which 
is outlined by the reddening, appears to be closely concentrated in its spiral arms. In Figure 
2 we present B and H greyscale images of the two galaxies. The B image shows that the 
foreground dust lane causes significant extinction of the background galaxy light, whereas 
in the H image the extinction is minimal, as expected. 



3. THE EXTINCTION CURVE 

3.1. Measuring the Extinction in Each Band 

To calculate the extinction through the disk of the foreground galaxy we use the 
method pioneered by [White fc Keel (1992)| . We call the surface brightness of the foreground 
spiral F, and that of the background spiral B. The observed surface brightness of the 
overlapping region will then be 

F + Be-^ (1) 

where r is the optical depth through the foreground spiral. The steps we follow to measure 
r for each band are as follows: a) we estimate B for regions where the two galaxies overlap, 
based on B in non-overlapping regions (more detail given in sec. 3.2); this estimate requires 
that the galaxy is symmetric; b) we estimate F by following a similar procedure as in step 
(a); c) we measure the surface brightness in the overlapping regions; d) we use eq. (1) to 
determine r as a function of position in the overlapping region of the foreground galaxy. 
More specifically, we measure the extinction Ax (where Ax = I.OSGta) in the spiral arm and 
inter arm regions. 



3.2. Spiral Arm and Interarm Regions 

We measure the surface brightness in square regions approximately 3."5 wide 
(1" = 170pc for a distance to the galaxy of 35Mpc), taking the mean pixel value to be 
our surface brightness estimate, and the uncertainty in the mean to be the error in our 
measurement. We do this for five different regions located along the foreground galaxy's 
outer spiral arm, where it occults the background galaxy. These five regions, which are 
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located at a radius ~ 1/3 from the center of the foreground galaxy, are marked in Figure 
2. We then estimate the background galaxy's true surface brightness by measuring it in 
five new regions, each of which is located at an opposite location (symmetrically reflected 
about the galaxy nucleus) from its parent region. Finally, we follow a similar procedure 
for the foreground galaxy to estimate its contribution to the total light seen in each of the 
overlapping regions. Due to the absence of a good counterpart for the foreground occulting 
arm, we estimate its average surface brightness from its own non-overlapping part. In this 
way, we obtain five estimates of the extinction through the foreground galaxy's spiral arm. 
Our final extinction estimate is the average of these five measurements. The errors are 
estimated by propagating the errors in each step of the above process. 

To measure the extinction in the interarm region, we repeat the same procedure, this 
time placing our first five squares in the interarm region just inside the foreground galaxy's 
outer spiral arm. These boxes, which are located at a radius ~ l.'l from the galaxy's center, 
are also marked in Figure 2. For comparison with the sample of occulting galaxies presented 
by [White et al. (1996)| , the B = 25mag/arcsec2 isophotal radius of NGC 2207 is Rg ~ 1.7. 



In these units, the regions that we examine are -RZ-Rfs = 0.76 for the spiral arm region and 
R/R2^ = 0.65 for the interarm region. 

Although our extinction measurements are independent of absolute calibration, we can 
check the quality of our galaxy color estimates by comparing them to other galaxy colors. 
The estimated colors of the foreground and background galaxies in the overlapping arm and 
interarm regions are listed in Table 1. Comparing these to galaxy colors found by [Frogel 



1988), we are reassured that they are reasonable colors for Sbc and Sc spiral galaxies. 



A good discussion of errors is given by [White et al. (199B}| . The greatest source of error 



in our extinction measurements is possible asymmetries in the background galaxy. Since 
the foreground galaxy is substantially fainter than the background galaxy, asymmetries in 
the foreground galaxy will not greatly affect our estimates of Ax. Small errors may also be 
introduced by imperfections in the flat-fielding and sky-subtraction of the images; however, 
these errors will not be as large as those caused by asymmetries in the background spiral 
galaxy. We estimate the average error in Ax that could be caused by such asymmetries 
by examining the background galaxy's deviations from symmetry in its non-overlapping 
regions. We find that the error in the measured extinction ranges from AAx ~ 0.08, in B, 
to AAx ~ 0.03, in K. The errors are smaller in the infrared bands because the background 
galaxy is much smoother (and, therefore, more symmetric) in the infrared than it is in 
the visible. The photometric errors, on the other hand, range from AAx, ~ 0.06 in B, to 
AAx ~ 0.03, in K. We sum these two errors in quadrature to obtain total error estimates 
for the extinctions we measure. 
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4. RESULTS 



Our measurements for the B VR J HK extinction through NGC 2207 in its spiral arm and 
interarm regions are hsted, along with their uncertainties, in Table 2; the total extinction 
errors that we find range from O.lmag, in B, to 0.04mag, in K. Like [White fc Keel (1992 



we find that the extinction is mainly concentrated in spiral arms (where Ay ~ 1.0), while 
the interarm regions (where Ay ^ 0.5) are mostly transparent. NGC 2207 has values for 
Ab at R/R25 similar to one of the dustier galaxies studied by [White et al. (199^ 



Also listed in Table 2 is the ratio F/B oi the foreground galaxy surface brightness 
to that of the background galaxy. This ratio indicates that the foreground galaxy's 
contribution to the total light in the regions examined is small. Our extinctions will, 
therefore, not be greatly affected by errors in our estimate of F. 

In Figure 3 we plot the extinction, normalized to the extinction in the J band, as a 
function of band wavelength. We are surprised to find that the extinction curve for the 
spiral arm region seems "fiatter" ("greyer") than expected from a standard Ry = 3.1 dust 
model, where Ry = Ay/E{B — V). In other words, the visual extinction is found to be less 
than expected compared to that in the near-infrared. The extinction curve appears to agree 
with an Ry = 5.0 dust model. This effect is even more pronounced in the interarm region 
extinction curve, where even an Ry = 5.0 dust model overestimates the visual extinction. 



For the bulk of the ISM in the Milky Way, the observed extinction law is Ry = 3.1; 
values of Ry = 5.0 are found in molecular clouds and HII regions ( [Savage fc Mat his 1979[ ; 



Uardelli et al. 1989| ). Deviations from this law, such as the one that we find, have been 



detected in various systems. One such system is the Sombrero galaxy, an edge-on system 
studied by [Emsellem (1995)] where stars and dust are mixed. Although [Emsellem (1995)] 



also observed an Ry = 5.0 extinction law, the geometry of the Sombrero galaxy is extremely 
different from that of our overlapping galaxies. Also, Bruzual et al., 1988 modeled the 



effects of dust mixed with stars and found that some deviations from an Ry = 3.1 law 
are expected due to scattering of light. [White et al. (1996][ calculated Ab/Aj for their 



overlapping galaxy candidates and found that some of them seem to have a flatter extinction 
curve than the Milky Way. Our more detailed study of the extinction curve conflrms this 
finding. 

Two effects may be responsible for the high Ry for the extinction curve that we 
observe. Scattering of light from the background galaxy into the line of sight can cause 
us to overestimate the amount of light successfully passing through the foreground disk 
and, thus, underestimate the extinction. This effect will be greater in the visual bands 
where scattering is more efficient, so a fiattening of the extinction curve is expected. An 
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unresolved patchy dust distribution also could be responsible for the high value of Ry- 



4.1. Scattering of Background Galaxy Light into the Line of Sight 



We can estimate the maximum effects of scattering by adopting dust grain properties 
and considering the geometry of our galaxies. We assume that the geometry consists of a 
simple dust screen (NGC 2207), perpendicular to our line of sight, that partially occults 
the background face-on galaxy (IC 2163). The scattering effects will be maximum if we 
consider the background galaxy to be a point source having a small scattering angle, 6. In 
this case, the maximum surface brightness in each band that could be scattered into the 
line of sight is approximated by 

Sma. = - (4.85x10-^)2 (2) 



Ivr \ a J 

where r is the optical depth, u is the dust albedo, Btot is the total brightness of the 
background galaxy, D is the distance to the background galaxy, a is the distance between 
the two galaxies, and last numerical factor converts S^ax into units of flux/arcsec^. Here we 
have assumed the Henyey-Greenstein phase function H{6) which, for ^ = 0, is 

H{0) = {l-g')il + g'-2g)-'/' (3) 

where g is the phase function asymmetry (we adopt = because this maximizes Smax)- 
Using N-body simulations and a distance to the pair of 35Mpc, p. Elmegreen et al.| 



1995) found that a = 40Kpc. So for our system 



Smax = 1.43 X IQ-^TLuBtotJ^^^^. (4) 

There is some uncertainty as to the values of uj and g for the near-infrared (especially at 
K). Many studies have suggested that the dust albedo and phase function asymmetry in 
the near-infrared are much higher than the predicted values of [Draine fc Lee (1984 )| (e.g.. 



Johnson et al. 199q ; [Pendleton et al. 1990| ; [Sellgren et al. 1992| ; [Witt et al. 1994t 



et al. l9^ ; Lehtinen fc IVlattila 199(:j ). In order to ensure that we are calculating the 



maximum effect of scattering, we use the high values for uj and g chosen by [Witt fc Gordon 



1996 )| . These values are cu = 0.6 and g = 0.6, in all bands. 



Table 3 lists the values we use to estimate Smax, the percent contribution of scattered 
light to the total light we see, and the maximum effect of scattering on our extinction 
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measurements. Our values for r may be found in Table 2 (where tx = 0.921 Ax). We 
estimate Bfot by integrating all the light from the right half of the background galaxy (the 
half which is not occulted by the foreground galaxy) and multiplying by two. Table 3 
lists our values for Btot and Smax in all observed bands. Table 3 also lists the value of the 
ratio Smax/{F + Be~'^) which indicates what fraction of the light we see in overlapping 
regions could be scattered light. Finally, Table 3 lists A Ax/ Ax, the maximum errors in our 
extinction estimates that could result from scattering of light. 

In all bands, AAx/Ax is less than 0.033. For purpose of comparison, AAx/Ax due 
to photometric and asymmetry uncertainties is greater than 0.073 in all bands, i.e., the 
combined photometric and asymmetry uncertainties (listed in Table 2) are much greater 
than the effects of scattering. Clearly, scattering cannot account for the high Ry value that 
we observe. Moreover, the effect of scattering on the interarm extinction curve is smaller 
than that on the spiral arm extinction curve, thus, making it even harder to explain their 
significant difference as due to scattering. 



4.2. Patchy Dust Distribution in NGC 2207 



Since scattering effects cannot explain the deviations from the Galactic extinction curve 
that we observe, we examine the possibility that this deviation is caused by unresolved 
clumping of the dust. It is reasonable to expect a patchy dust distribution in NGC 2207. 
The dust distribution in our own galaxy is patchy on scales small enough that the patches 
would be unresolved if seen from a distance of 35Mpc ([Savage fc Mathis 1979|) . Moreover, 



Rix fc Rieke (1993)| have mapped the dust distribution in M31 and have shown that it is 



quite patchy in appearance. [White et al. (1996)] suggest that dumpiness of dust is probably 



the cause for their overlapping galaxies' flat extinction curves. [Witt fc Gordon (1996)[ have 
studied the radiation transfer in a two-phase clumpy dust medium and have also suggested 
that a clumpy medium tends to flatten the effective extinction curve. 

Is a patchy dust distribution in NGC 2207 responsible for flattening our derived 
extinction curve? To investigate this we explore a simple two-dimensional, two-component 
dust model. We assume that the dust in the foreground galaxy consists of a high density 
component and a low density component. The high density dust component has a J-band 
optical depth ri, and an area filling factor /. The low density component has a J-band 
optical depth T2 and an area filling factor 1 — /. The effective transmission through the 
patchy dust distribution will then be 



+ (1 _ /)e-^ (5) 
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where Tgff is the effective optical depth that we measure. We assume that each dust 
component follows an Ry = 3.1 Galactic extinction law, and study the extinction curve for 
Teff- We find that any combination of values for our free parameters Ti, T2, and / tends 
to flatten the effective extinction curve. In particular, the slope of the effective extinction 
curve appears to be most sensitive to differences in ti/t2 and /; the curve is flatter for 
larger values of ti/t2 and for values of / which are closer to 50%. 

For purposes of illustration, we choose two sets of parameters, one for the spiral arm 
and one for the interarm, that reproduce the effective extinctions that we measure (and 
thus, the extinction curve that we observe.) These parameters are: ti = 2.0, T2 = 0.22 and 
/ = 0.12 for the spiral arm region, and Ti = 2.0, T2 = 0.08 and / = 0.15 for the interarm 
region. These models give effective J-band extinctions of A^ff = 0.35 and ^e// = 0.22, 
respectively, which match the measured extinctions in Table 2. We plot these models 
in Figure 4 to illustrate that the two-component dust model is capable of reproducing 
our extinction curves. We are cautious about assigning any physical significance to these 
specific models however, because their uniqueness is not established. 

In order to fully explore the parameter space of our two-component model, we compute 
for our measured extinctions, keeping / fixed and varying ri and T2. In Figure 5, we show 
contours of constant reduced-x^ for four different values of the filling factor /. Contours 
are shown for both the spiral arm (bold) and the interarm (light) extinctions. We note that 
T2 is reasonably tightly constrained and does not show a significant dependence on ti or 
/. Ti, on the other hand, is not that well constrained. The "golf club" contours in Figure 
5 reveal the main difference between the spiral arm and interarm regions: the low density 
dust component, T2, is smaller in the interarm region than it is in the spiral arm region. In 
other words, the observed arm-interarm difference can be explained by a model containing 
high density dust clumps embedded in a low density interclump medium which is denser 
in the galaxy's spiral arms, and less dense in its interarm regions. An example of this 
model would be molecular clouds (corresponding to ri) embedded in a lower density dusty 
environment (corresponding to T2). One could then speculate that, as the ISM moves from 
an interarm into a spiral arm region, the low density interclump regions are compressed {t2 
increases), whereas the molecular clouds, being bound systems, are not much affected (ri 
stays approximately the same). According to the two-component patchy dust model, this 
scenario may be able to explain the extinction laws that we observe. 

In addition to our two-component patchy dust model, we also explore the clumpy dust 
model of |Witt fc Gordon (1996)|. This is a three-dimensional two-phase model. The model 



consists of a spherical volume divided into bins, each of which is randomly assigned 
to either a high-density state or a low-density state. Photons originate from a centrally 
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located, isotropically radiating point source and their transfer through the clumpy medium 
is studied. Parameters of this model are N, the volume filling factor //, the density 
ratio between the low-density and high-density phases ^2/^1, and the optical depth of the 
equivalent homogeneous dust distribution th- The model is described in great detail by 
[Witt fc Gordon (1996)| . The main advantage of this model over our two-component model 



is that it is three-dimensional, thus allowing for each line of sight to have a variety of optical 
depths (depending on the number of high-density bins that are encountered). Although its 
geometry differs from that of our overlapping galaxy pair (i.e., it contains a point source 
of light embedded in a spherical dust distribution, rather than an extended light source 
behind a flat dustscreen), the [Witt fc Gordon (1996)| model is, in this case, essentially 



indistinguishable from a similar model of planar geometry since scattering is negligible. 

We fit this clumpy dust model to our measured visual and infrared extinctions by 
minimizing x^- We perform the fits for the case = 10, and our free parameters are 
^2/^1 and the J-band value of th- The fits are relatively crude because we do not explore 
the full parameter space of the model; nevertheless, they suffice to demonstrate that the 
Witt fc Gordon (1996)| clumpy dust model is capable of producing our observed extinction 



curves. In Figure 6 we present the best fits for the spiral arm and interarm extinction 
curves. The best-fit values of our parameters are /C2//C1 = 0.032, th = 0.58 for the spiral 
arm region, and k2/ki = 0.0032, th = 0.50 for the interarm region. These models give 
effective J-band extinctions of A^ff = 0.35 and A^ff = 0.24, respectively, which match 
the measured extinctions in Table 2. These results are in agreement with those of our 
two-component model, as they attribute the observed difference between the spiral arm and 
interarm extinction curves to the interarm region having a higher degree of dust patchiness 
(density ratio between high-density and low-density phases). 

The fits are not perfect (especially in the spiral arm region), but this is not unexpected 
because the full parameter space of the model has not been explored. Future clumpy 
dust models, suggested by [Witt fc Gordon (1996)[ , that examine the case of an extended 
star distribution will be particularly well suited to explain the behavior of our observed 
extinction curves. 



4.3. The Column Depth of Hydrogen in NGC 2207 



Besides flattening the extinction law, an unresolved patchy dust distribution will cause 
us to underestimate the total column depth of hydrogen in a galaxy if we use the visual 
extinction for this purpose. To investigate this we use the two-component model discussed 
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in sec 4.2. According to this model, the ISM is composed of two components, a high density 
one and a low density one, having J-band optical depths of ti and r2, respectively. The 
effective extinction through this ISM is given by equation (5), where / is the area filling 
factor of high density clumps. The optical depth that corresponds to the total column 
depth of hydrogen, however, is given by 

TH = fn + (1 - /)r2 (6) 

where th is directly proportional to the column depth of hydrogen, Nh- The relation 
between Nh and the F-band extinction Ay is given by ([Savage fc Mathis 1979| ) as 
Nh ~ 1.87 X lO^^atom cm-'^Av- 

Adopting the values of ti, T2, and / shown in Figure 4, and assuming an Ry = 3.1 law, 
we find that the total column depth of hydrogen in NGC 2207 is: 3.1 x 10^^ atoms cm~^ in 
the spiral arm region, and 2.7 x lO^^atoms cm~^ in the interarm region. If we had used the 
observed F-band extinctions (listed in Table 2) to estimate Nh we would have found a total 
column depth equal to 1.8 x lO^^atoms cm~^ and 1.0 x lO^^atoms cm~^ in the spiral arm 
and interarm regions, respectively. If, on the other hand, we had used the observed K-hand 
extinctions to estimate Nh, we would have found a total hydrogen column depth equal to 
3.8 X lO^^atoms cm~^ and 2.4 x lO^^atoms cm~^ in the spiral arm and interarm regions, 
respectively. In other words, when using the F-band extinction, we underestimate the total 
column depth of hydrogen by a factor of ~ 2 solely as a result of an unresolved clumpy dust 
distribution. In contrast, the iT-band extinction provides column depth estimates that are 
much more reliable. 

The above example clearly demonstrates that it is much better to use the infrared 
extinction in order to estimate the total column depth of hydrogen in a galaxy. The reason 
for this is that most of the hydrogen is in dense components which are more easily seen in 
the infrared than in the visual. 



5. SUMMARY & DISCUSSION 



The conclusions of this paper are as follows: 



1. The dust in NGC 2207 is mainly concentrated in its spiral arms, leaving its interarm 
regions mostly transparent. This confirms the results of [White et al. (1996)| for other 
overlapping galaxy pairs. 
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2. The extinction curve in the spiral arm region is flatter than expected, as it resembles a 
Milky Way Ry = 5.0 dust model. This conflrms the suspicions of [White et al. (1996 



3. The extinction curve in the interarm region is signiflcantly flatter than that of the 
spiral arm region. This is, perhaps, the most interesting of our results and could be 
a valuable clue to the nature of the dust distribution in different regions of galaxies. 
It would be very interesting to know whether this effect is observed in other spiral 
galaxies. Similar studies should be performed on other overlapping galaxy pairs in 
order to answer this question. 

4. Scattering effects are negligible in this overlapping galaxy pair. 

5. An unresolved patchy dust distribution in NGC 2207 is capable of producing the 
extinction curves that we observe. This conflrms the suspicions of [Witt fc Gordon 



(1996)1 and [White et al. (1996)| . Since we are averaging over boxes which are ~ 600pc 



on a side, we expect clumping of dust to be signiflcant on these size scales. Studies 
that probe the ISM on much smaller scales might be less susceptible to patchy dust 
effects. 

Fits of two patchy dust models to the data suggest that the arm-interarm difference 
in the observed extinction law can be explained if the galaxy's interarm region has a 
higher degree of dust patchiness (density ratio between high-density and low-density 
phases) than its spiral arm region. 

It is clear that unresolved clumping of dust will cause one to underestimate the 
average column depth of gas in a galaxy if based solely on the visual extinction. It is 
much better to use the infrared extinction for this purpose. 



It would be interesting to compare our derived extinction laws for NGC 2207 with 

[Rix & Rieke (1993)[ use the 



those that might be obtained from the galaxy's colors (e.g. 
colors of M51 to determine its optical depth). Also, future efforts should be made to obtain 
high resolution visual and infrared images of very distant overlapping galaxies. The method 
developed by [White et al. (1996][ , which makes use of overlapping galaxy pairs, is a valuable 
tool for probing the ISM in distant galaxies. It is, perhaps, one of the few ways to study 
the dust content in high redshift galaxies. 
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TABLE 1. Derived Intrinsic Galaxy Colors 



color 


NGC 2207 


IC 2163 




arm 


interarm 


arm 


interarm 


B-V 


0.52 


0.71 


0.80 


0.66 


V-R 


0.45 


0.52 


0.63 


0.56 


R-J 


1.25 


1.30 


1.80 


1.64 


J-H 


0.77 


0.67 


0.66 


0.63 


H-K 


0.07 


0.29 


0.30 


0.27 



Note — Estimated errors range from 0.04 mag in Kto 0.10 mag in B. 
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TABLE 2. The Extinction through NGC 2207 





Spiral arm 


Interarm 




Filter 


Ax 




F/B 


B 


1.32 ±0.10 


0.60 ±0.10 


0.27 


V 


0.96 ±0.10 


0.52 ±0.10 


0.21 


R 


0.80 ±0.09 


0.44 ±0.09 


0.17 


J 


0.35 ±0.06 


0.24 ±0.07 


0.10 


H 


0.27 ±0.05 


0.19 ±0.06 


0.11 


K 


0.22 ±0.04 


0.14 ±0.06 


0.09 



Notes — ^4;^ is the extinction measmcd at a wavelength A. F/B is the ratio of the foreground 
galaxy surface brightness to that of the background galaxy. Uncertainties shown are 
photometric uncertainties and uncertainties due to asymmetry in the background galaxy 
added in quadrature. 
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TABLE 3. 


. Scattering 


; Estimates 




Filter 


Btot 


Smax 


Smax 

F+Be—^ 


AAx 
Ax 




(mag) 


(mag/D") 






B 


13.03 


25.49 


4.0% 


3.3% 


V 


12.28 


25.07 


2.4% 


2.7% 


R 


11.70 


24.71 


1.7% 


2.4% 


J 


10.00 


23.90 


0.5% 


1.7% 


H 


9.35 


23.52 


0.4% 


1.6% 


K 


9.10 


23.51 


0.3% 


1.5% 



Notes — Btot is the total brightness of the background galaxy; S,nax is the maximum surface 
brightness that could be scattered into the line of sight; F and B are the surface brightnesses 
of the foreground and background galaxies, respectively; r is the extinction optical depth 
through the disk of the foreground galaxy; A\ is the extinction of light (as a function of 
wavelength) through the disk of the foreground galaxy; AAx is the maximum change in our 
extinction estimates that could result from scattering of light. 
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FIG. 
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late 



Fig. 1.— Color image of NGC 2207 (right) and IC 2163 (left) constructed from HR and 
i?-band images. The image has dimensions (6.'9 x 8.'6). Note that IC 2163 can be easily seen 
through NGC 2207. Also note the reddening of background galaxy light in the foreground 
galaxy's spiral arm region. 




60" 



Fig. 2.— B and i/ images of NGC 2207 (right) and IC 2163 (left). In these images North 
is up and East is to the left. The square regions in which we measure the extinction are 
marked on the B image (there are five in the spiral arm and five in the interarm region). 
The regions are approximately 3."5 wide (1" = 170pc for a distance to the galaxy of 35Mpc). 
The angular scale of the images is given at the bottom. Note that the B image shows that 
the foreground dust lane in NGC 2207 causes significant extinction of the background galaxy 
light, whereas in the H image the extinction is minimal. 
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Extinction Curves for NGC 2207 
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Fig. 3. — Extinction vs. wavelength for the dust content in NGC 2207's spiral arm (top) 
and interarm (bottom) regions. The points represent our extinction measurements Ax, as 
a function of wavelength A, at BVRJHK, normalized to Aj, the extinction in J, along with 
their la uncertainties (a discussion on errors can be found in sec. 3.2). The solid and 
dashed lines represent Ry — 3.1 and Ry — 5.0 Galactic extinction laws respectively (where 
Ry — Av/E{B — V)). Note that the extinction laws in both the spiral arm and interarm 
regions of NGC 2207 are flatter ("greyer") than a standard R — 3.1 Galactic extinction law. 
Also note that the interarm extinction law is substantially flatter ( "greyer" ) than the spiral 
arm extinction law. 
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Two Component Patchy Dust Model 
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Fig. 4. — Effective extinction laws for a two-component dust model. The two panels show 
the effective extinction curves tailored to fit the spiral arm (top) and interarm (bottom) 
extinction curves. The points represent our extinction measurements Ax, as a function 
of wavelength A, at BVRJHK, normalized to Aj, the extinction in J, along with their 
la uncertainties. The solid lines represent an Ry = 3.1 Galactic extinction law (where 
Rv = Ay/E{B — V)). The dashed lines represent the effective extinction curves which 
result from our patchy dust model. The values of ti, T2 (the J-band optical depths of high 
and low density dust components, respectively), and / (the area filling factor of high density 
clumps) corresponding to these curves are shown. These parameters are defined in sec. 
4.2. Note that the two-component model does an excellent job in reproducing our observed 
extinction curves (this figure is meant to illustrate how a patchy dust distribution is capable 
of explaining our observed extinction curves. The specific models shown above are not unique 
and should not be assigned much physical meaning.) 
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Fig. 5. — Contours of constant for our two-component dust model. The four panels show 
contours of constant reduced-x^ as a functionof the paramet ers Ti and T2 (the J-band optical 
depths of high and low density dust components, respectively). The contours represent fits 
of this model to our observed spiral arm and interarm extinctions for area filling factors of: 
/ = 0.12 (top left), / = 0.15 (top right), / = 0.20 (bottom left), and / = 0.25 (bottom 
right). Bold and light contours represent fits to the spiral arm and interarm extinction 
curves, respectively. The contours shown are for values of the reduced-x^equal to 3, 2, 1, 
and 0.5 (going from outer to inner contours). Where only three contours are shown, these 
are for values of the reduced-x^ equal to 3, 2, and 1. We note that T2 is reasonably tightly 
constrained and does not show a great dependen ce on ti or /. ti, on the other hand, is not 
tightly constrained. These "golf club" contours reveal the main difference between the spiral 
arm and interarm regions: the low density dust optical depth, T2, is smaller in the interarm 
region than it is in the spiral arm region. 
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Witt & Gordon (1996) Patchy Dust Model 
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Fig. 6. — Effective extinction laws for the Witt & Gordon (1996) clumpy dust model. The 
two panels show the best-fit effective extinction curves for the spiral arm (top) and interarm 
(bottom) extinction curves. The points represent our extinction measurements Ax, as a 
function of wavelength A, at BVRJHK, normalized to Aj, the extinction in J, along with 
their la uncertainties. The solid lines represent an Ry = 3.1 Galactic extinction law (where 
Rv = Ay/E{B — V)). The dashed lines represent the best-fit effective extinction curves 
which result from the Witt & Gordon (1996) clumpy dust model. The best-fit values of 
k2/ki (the density ratio between low density and high density phases) and th (the J-band 
optical depth of the equivalent homogeneous dust distribution) corresponding to these curves 
are shown. Note that this model suggests that the observed difference between the spiral 
arm and interarm extinction curves is caused by the interarm region having a higher density 
contrast between high and low density clouds {ki/k2) than the spiral arm region. 



